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A free field representation for the type / vertex operators and the corner transfer matrices of the eight- vertex model 
is proposed. The construction uses the vertex-face correspondence, which makes it possible to express correlation 
functions of the eight-vertex model in terms of correlation functions of the SOS model with a nonlocal insertion. 
CN I This new nonlocal insertion admits of a free field representation in terms of Lukyanov's screening operator. The 

spectrum of the corner transfer matrix and the Baxter-Kelland formula for the average staggered polarization have 
been reproduced. 
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1. Introduction 

The eight-vertex model 1,2 ' 3 is one of the most famous examples of exactly solvable models of statistical 
mechanics. Its partition function was found in 1971 by Baxter. 4 Later expressions for the spontaneous 
polarization and magnetization were proposed. 5,6 Nevertheless, very few results on correlation functions 
have been obtained up to now. 3,7,8 
Qh. For the last years, a new approach to calculating correlation functions of integrable (exactly solvable) lat- 

tice and continuous models has been developed. 9-13 This approach, known as the vertex operator approach, 
has shown its efficiency in calculation and analysis of correlation functions and form factors for different inte- 
grable models, such as the six-vertex model, 9-11 the restricted solid-on-solid (RSOS) models 14-16 and their 
An generalizations, 17,18 the fusion RSOS models, 19,20 the Ising model, 14,21 the sine-Gordon model. 12 ' 13 
The most important constituent of this approach that allows one to perform calculations is the free field 
representation of the vertex operators. 

However, the free field representation was not elaborated for the eight-vertex model. Formally, the 
problem lies in the presence of the d matrix element in the R matrix of this model which breaks down the 
charge conservation, which is necessary for the free field representation. In principle, the same reason gives 
no way of applying directly the Bcthe ansatz method to this model. 3 Baxter's idea to overcome this difficulty 
for the Bethe ansatz was to twist somehow the model to restore the charge conservation law without change 
of eigenvalues of the transfer matrix. 22 Such 'twisted' model was the solid-on-solid (SOS) model, which is a 
face type model. It is related to the eight-vertex model by the so-called vertex-face correspondence. 

It is an old idea to use this vertex-face correspondence to relate correlation functions of the eight- 
vertex model with ones of the SOS model, 23 and then to apply the known free field representation for the 
vertex operators of the SOS model. The obstacle is non-local character of the vertex-face correspondence 
for correlation functions. In this paper we analyze this relation and propose a free field approach for the 
correlation functions of the eight-vertex model. As a check we find the spectrum of the corner transfer 
matrix of the eight-vertex model, and reproduce the Baxter-Kelland formula 5 for the staggered spontaneous 
polarization in the antiferroelectric phase. 

The paper is organized as follows. In Sec. 2 basic definitions and notations are given. In Sec. 3 we 
give a review of the vertex operator approach of both the eight-vertex and the SOS models. In Sec. 4 we 
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Fig. 1. Eight-vertex model: (a) definition of the weight matrix; (b) two degenerate ground states. 



use the vertex-face correspondence to find a relation between correlation function of the eight-vertex model 
and the SOS model. The construction contains a new object in addition to corner transfer matrices and 
vertex operators. This object (operator A in our designations) relates the corner transfer matrices of the 
eight-vertex model to those of the SOS model. In Sec. 5 we describe a free field representation and calculate 
the staggered polarization by use of the free field representation. In Sec. 6 we discuss main conjectures used 
in the paper and possible directions of future studies. 
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2. Basic Definitions and Vertex-Face Correspondence 

2.1. Eight-Vertex Model 

The eight-vertex model is defined as follows. 3 The fluctuating variables (polarizations) e = ±1 = ± are 
placed at links of a square lattice. Interaction is associated with vertices of the lattice. A local weight R%\%* 
is associated to each configuration of polarizations e\, £2, £3, £4 around a vertex (Fig. la). A weight of a 
lattice configuration is the product of local weights. 

It is convenient to attach a constant complex 'spectral parameter' to each line of the lattice, and to 
consider a weight as a function of the difference of the spectral parameters of the two lines that intersect at 
the vertex. The weights as functions of this difference are given by 3 



(2.1a) 



with 9i(u; r) is the zth theta function with the basic periods 1 and t (Imr > 0). Here e and r are parameters. 
It is often convenient to use another set of parameters 

x = e- f -, p = x 2r , z = x 2u . (2.2) 

We shall use both 'additive' parameters u, e, r and 'multiplicative' parameters z, x, p on equal grounds. 
The common factor k(u)Rq(u) is chosen so that the partition function per site is equal to one: 

«(«) = z-^x 1 -i(^;^)- 2 (x 4 -^)^(x 2 z- 1 ;x 2r )- 1 ( 2 ; 2 '-- 2 z;x 2 '-) M 1 , (2.1b) 
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ni,...,njv=0 

Splitting the common factor into two parts is convenient for later use. 
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The matrix R(u) satisfies the Yang-Baxter equation, the crossing and unitarity properties. 3,8 This 
supplies commutativity of transfer matrices with different spectral parameters u (but coincident e and r) 
and integrability of the model. 

For real parameters e, r, u in the region 



e > 0, r > 1, -1 < u < 1 or a + |6| + |d| < c 



(2.3) 



the model is in the antiferroelectric phase. In the 'low temperature' limit e — > oo (or a, |6|, |d| <C c) the 
system falls in one of two ground states (Fig. lb). These two ground states correspond to the broken Z 2 
symmetry. At finite temperatures the symmetry remains broken in the whole region (2.3). The second order 
phase transition takes place at e = 0. The weight cLa in our definition is negative but it enters into any 
configuration in an even power in the antiferroelectric phase and its sign is inessential. 

The aim of the present paper is to find a procedure for calculating correlation functions in the eight- 
vertex model in the thermodynamic limit. In our normalization the partition function per site is equal to 
1. So we can define the correlation function directly on the infinite lattice. Let us label the vertices of the 
lattice of finite size by two integer Cartesian coordinates (x,y). The coordinates of polarizations e(x,y) arc 
the coordinates of the centers of links (x £ Z + i, y £ Z for vertical links and x £ Z, y £ Z + | for horizontal 
ones). The probability that polarizations at definite links (xk,yk), k = 1, . . . ,N, take fixed values e\. is given 

by 

p« „ = V n R ( u * ^ {x ' y+i)£{x+hy) 

x,yeZ 



E 



v v) 



(2.4) 



Here the sum is taken over all configurations such that £ Xk ,y k — £fc and e x . y only differ from the ground state 
configuration value (— ) x ~y~ 2+* at a finite number of links. Note, that we consider an inhomogeneous lattice, 



where each line brings its own spectral parameter u x of v y , so that u 



const and v v 



const for 



large enough x and y respectively. Later we concentrate our attention on the particular case of N parallel 
neighboring links (xk = 5, J/fc = k, see Fig. 5a), but our consideration is valid in general case. 



2.2. SOS model 

Let us turn to the SOS model. Here one associates an integer variable ('height') n to each site of a square 
lattice, so that the heights m and n 2 at adjacent sites satisfy the admissibility condition 



K - n 2 \ = 1. 



(2.5) 



Interaction is described by weights W 
attached to lines of the dual lattice (Fig. 2): 
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Note that the weights W contain [n] in the denominators. So the construction demands some regularization. 
Namely, we may think that the heights n belong to the set Z + S with some real S. Such shift does not 
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Fig. 2. SOS model: (a) definition of weights; (ft) an infinite number of degenerate ground states labeled by an integer 
m and by s = ±1 satisfying (k — l)r < m, m + s < kr for some integer k. 



break down the admissibility condition (2.5). After calculation of all physical quantities one may take the 
limit 8 — > 0. The consideration below is valid literally for the regularized case, and we shall imply this 
regularization everywhere. Moreover, the final results do not depend on 6 at all. 

The weights (2.6) satisfy a face version of the Yang-Baxter equation, the crossing and unitarity proper- 
ties. 22 ' 14 Note that the elements (2.6a) are analogues of the element a(u) of the eight-vertex model, (2.6b) 
are analogues of c(u), and (2.6c) of b(u). Namely, we can associate a sign variable e to each link of the dual 
lattice intersecting a link of the direct lattice between n + e and n (n + e is to the left or upper than n) . It 
is impossible to imagine any analog of d(u) in this picture. The charge defined as the sum of the variables e 
along an infinite line consisting of links of the direct lattice is conserved by definition. 

We are interested in the so called regime III: 

e>0, r>l, 0<u<l. (2.7) 
There is an infinite number of degenerate, but inequivalent, ground states in this regime shown in Fig. 2b. 

2.3. Vertex-Face Correspondence 
Consider the functions (Fig. 3a) 

, x„' e " _m ( n f(n'-n)u + n' .2ir\ n f(n'-n)u + n' .2tt\\ , . 



with m being an arbitrary integer (in the regularized version m € Z + 5), e = ± being a polarization variable, 
n and n' 
(Fig. 4a) 



n and n' being an admissible pair of heights (i. e. \n' — n\ = 1). These functions enter into the identity 22 
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which is referred to as vertex- face correspondence. The functions t e (u)™ are called intertwining vectors (the 
subscript e is often regarded as a vector index). Note that the identity holds for arbitrary value of uq. The 
vertex-face correspondence is the basic relation for solving the eight-vertex model by Bethe ansatz. Below 
we apply it to find the free field representation. The intertwining vectors may be pulled through the lattice 
by use of the vertex-face correspondence." 

It is convenient to introduce a 'conjugate' and a 'primed' intertwining vectors (Fig. 3b) 

E'e'W'M");- = <W, £ WXM"" = <W. (2.10) 



a The standard vertex-face correspondence, as it is described in Ref. 22, relates the SOS model in the regime III to the 
eight-vertex model in the disordered region. But we can relate the disordered and ferroelectric region by use of the Baxter 
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Fig. 3. Intertwining vectors: (a) graphic notations; (6) definition of the conjugate and 'primed' vectors. 
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Fig. 4. Vertex-face correspondence: (a) by usual intertwining vectors; (ft) by conjugate intertwining vectors. 
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The conjugate intertwining vectors enter into another form of the vertex- face correspondence (Fig. 4b): 
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3. Vertex Operator Approach 

In the vertex operator approach (see Ref . 1 1 and references therein) the computation of correlation functions 
is reduced to the algebraic problem of finding the traces of some operators. Namely, the lattice is partitioned 
into several parts as it is shown in Fig. 5. Then the partition functions of the parts with fixed variables 



duality transformation 3 

a{u) = ± (a D («) + b D (u) + c D («) + d D («)), 
b(u) = \{-a D (u) - b D (u) + c D (u) + c D (u)), 
c{u) = i (a D (u) - b D (u) + c D (u) - d D (u)), 
d(u) = i(-a D (u) + b D (u) + c D (u) - d D (u)), 

where ajj, bjj, cp, djj mean the weights in the disordered phase. This duality for correlation functions 24 gives 

e~n—m 

te(u)Z =-^=-(t£(u)» +st D (u)Z). 
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Fig. 5. Partition of the lattice into the corner transfer matrices A, B, C, D and vertex operators <E>. The lattice 
is only cut to stress the mode of partition, it must be considered as an integral thing. Action of the operators is 
anti-clockwise, (a) Eight-vertex model, (ft) SOS model. 

at their boundaries are treated as matrix elements of some operators (corner transfer matrices 3 and vertex 
operators 9 ) acting in the spaces associated with a half-line. 

3.1. Eight- Vertex Model 

Let us begin with the eight-vertex model in the antiferroelectric region. 8 Consider the subspace Hi in 
C 2 (£> C 2 (£> • • •, spanned on the vectors \s\) ® \e2) ® • • • such that Ek stabilizes to (— ) k+% for large k. The 
operators associated to the parts of the lattice, shown in Fig. 5a, are 

^-^(uJ.C^-^^u) : Hi -» Hi-i, B^{u),D^{u):Hi^Hi, 
^-^(u),^ 1 "^'^) : H % -» 
The product of four corner transfer matrices in the infinite lattice is an operator independent of u: 

pW = ^ consta; 4ff(,) , (3.1) 

where const is a c-number and if W i s so called corner Hamiltonian. 

The spectrum of the operator is equidistant. More precisely, it is determined by the generating 
function 

^WeT^/^^u); 1 . (3.2) 

The vertex operators <j>i 1_J '^ (u) satisfy the following relations which are consequences of the Yang- 
Baxter equation, unitarity and crossing symmetry: 



(3.3) 



^-^(uJpW = p( 1 -*)*( 1 -*'*)( ti _ 2), 



(3.4) 
(3.5) 



The conjugate vertex operator (J^ 1 ' 1 ^*(u) denned as 



$( 1 - i ' i ) / (u)B (i) (u)4 (i,1 ~ i) («) = B^i^A^-^iu)^' 1 -^*^), E $£ <,1_<) *(u)*i 1_i,<) (u) = 1 (3.6) 

s 

is given by 

^-i)*( u ) = ^- i \ u -l). (3.7) 

The probability of the configuration of polarizations on an inhomogeneous lattice as it is shown in Fig. 5a is 
given by 



Consider more general functions 

F^.. £N (u 1} ...,u N ) = Tr Wi (d^W) . . . ag-^tii) x 4 ^') (3.9) 
for even TV. They satisfy the equations 8 

O-rri 

41™ (• ..,«,• + —,...) = 41e„(- ■■,«*■■ ■), (3-lOa) 

^ ( i } .. £w («i + o, . . . , ujv + o) = FW.. £jv ( Ul , . . . , UJV ) (3.10b) 

53 - »i- -.,...(• . .,uj,u j+1 , . . .) 

£? ' +1 =F% +i£j ^...,u j+1 , Uj ,...), (3.10c) 

^ ( ;k... £jv K + 2,«2, ...,un)= Fg-% ei {u 2 , . . .,u N , Ul ), (3.10d) 
E *£L w , e ,- e («i> ■ • • « - 1) = ■ ■ ■ . E F -U u > «-!) = !■ (3-10e) 



These equations follow from Eqs. (3.4-3.5) and can be considered as defining equations for correlation 
functions. 



3.2. SOS model 

For the SOS model the construction repeats that for the eight-vertex model in general outline (see Fig. 5b). 
The corner transfer matrices are labelled by two integers m and n. The variable n denotes the height at the 
corner. The number m describes the condition at the infinity. Namely, for the space H m ,n with < in < r — 1 
is the subspace of C°° <g> C°° <£> • • • spanned of vectors |n ) ® \ni) <g> |n 2 ) <g> • • •, such that n = n, each pair 
(rife, rik+i) is an admissible pair of heights, and for k large enough n 2 k = m and n 2 k+i = m + 1 if n — mis 
even, n 2 fe = m + 1 and n 2 fc+i = m if n — m is odd. The situation of general m is more complicated, but 
we need not discuss it for our purposes. For simplicity we restrict the discussion to r > 2, but our final 
construction is valid for any r > 1. Parts of the lattice shown in Fig. 5b define again operators: 

-^■m,n(^), -^m,n(^), ^m,n(^0, -^m.n(^) • ^m.n * ^m,n, 

The product of the corner transfer matrices is again a u-independent operator: 

p m . n = D m . n (u)C m ,n(u)B m . n (u)A m . n (u) = COnSt [ft] X iHm -" . (3.11) 
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The operator H m>n also possesses an equidistant spectrum. We assume that the respective generating 
function is 

XmAl) = Tr « m ,„ q Hm - = q(™-('-l)«) 2 /4r(r-l) ((Z; (3^) 

The vertex operators 3>(m)^ (we usually omit indices m in operators from here on) obey the relations: 
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s n 





d>M£cD( Wl )*, (3.13) 



§{u) n nPn = Pn ,$>{u-2) n n , (3.14) 
£(n' - n)[n') ${u - 1)^$«' = (-1)-™. (3.15) 

n' 

The conjugate operator defined as 

&(utfB n (u)A n (u) = B n ,{u)A n ,{u)&{u)i, £ = 1 (3.16) 

is given by 

- (-)""">' - n)[n] *(« - ltf. (3.17) 
4. Vertex-Face Correspondence and Correlation Functions 

Now we are ready to explore the vertex-face correspondence for the vertex operator construction. 

Let us enclose a piece of lattice of the eight-vertex model with several fixed polarizations with the 
intertwining vectors from the upper and right sides and with the conjugate intertwining vectors from the 
lower and left sides (Fig. 6a). We can shrink such loop using the vertex-face correspondence and the definition 
of the conjugate intertwining vector (2.10) (Fig. 6b). But the tail in the right lower part of Fig. 6b cannot 
be cancelled, because the heights n and n' designated at the figure are generically not equal. 

Now let us establish correspondence between ground states (or conditions at infinity) for the SOS and 
eight-vertex models. Let m in the definition of the intertwining vectors (2.8) coincide with m characterizing 
the condition at the infinity. To do it, consider the 'low temperature limit' x — > (e — > oo). For simplicity 
we consider the case < m < r — 1. 

It can be checked (see Appendix A) that the relevant values of u in t e (u)™ and t*(u)™ for the free field 
representation are 

-2 < Rcu < 0. 
In the limit x — > we have for the intertwining vectors 

t+(«)n ~ ^e-^"'-"^"') 2 , 

t (it)™' ~ ( — )™~" l + 1 / £r (("'-")"+"') 2 -tr+e((n' -n)u+n') 

' " ^ ^ ' (4.1) 

t+(u)n ~ (n - n') J— e ^(( n -™')"+n) 2 + (l+™'-n)e« ; 

t*_{u)i ~ (-) n - m +\ n ' - n )^ e i- r ({™-n')u+n) 2 +er-e(n-u)^ ^ > q 

Consider now the quantities (Fig. 7a) 

a£-^ =t* e {u)2V±- i tM2X\- i - (4-2) 

It is easy to check that 

|a^ 0) | ~ 1, |a<_ 1>0) | ~ e 2eRe «, \a^\ ~ e 2 ^«, \a^\ ~ 1. 
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Fig. 6. Contraction of the loop of intertwining vectors in the lattice with insertions, (a) Initial loop. The t' vector 
in the lower write side is used to simplify the next figure. (6) 'Maximal' contraction. Dotted lines connect points 
with coincident heights. Summation in n and n' is implied, so they do not necessarily coincide, (c) Partition of the 
lattice after such contraction. 
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Fig. 7. Vertex-face correspondence of ground states: (a) basic weights; (ft) configurations of maximal weight at low 
temperatures (x — » 0) for m > 0. 



We see that 

|a^|»|aL 1 '° ) |, |aL°' 1) |»|aV°' 1) |. 

The configurations shown in Fig. 7b give the leading contribution into the partition function. The other two 
configurations are vanishing in the low-temperature limit. It means that a definite ground state configuration 
of the SOS model imposes a definite ground state configuration of the eight-vertex model. Hence the SOS 
condition at the infinity m, such that n — m E 2Z + i, imposes the condition at the infinity i for the 
eight- vertex model. 

So we can partition the lattice as it is shown in Fig. 6c. According to this partition we introduce 
the following objects: the corner transfer matrices of the SOS model, the vertex operators, and the 'tail 
operator'. 
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Let us begin with the corner transfer matrices of the SOS model 



p'W = Dm, n (u)C m , n (!i)Bm,n(«)Aii,n(M) = Const [n] X 
ne2Z+m+i ne2Z+m+i 



4H n 



(4.3) 



We assume that one have to keep fixed an arbitrary nonzero value of m to calculate correlation functions. 
This conjecture seems to be physically reasonable, because such condition is compatible with the vertex-face 
correspondence: we may take respective values of heights at the boundary of Fig. 6a, not disturbing the 
shrinking procedure. Besides, the correlation function in the eight-vertex model must be sensitive to the 
eight- vertex boundary condition labeled by i = n — m (mod 2), but insensitive to details of arrangement of 
the facilities that create these conditions, like a system of intertwining vectors and the SOS lattice somewhere 
at infinity. 

The vertex operators are given by 



<s>i i -^(u:u ) = f(u-u a ) t e (« -«„);:'*«', 

$( 1 - i - i )*(«;u ) = r 1 («-«o) t:(«-uo)s'**«'=^r M) («-i) 



(4.4) 



with the function 



/(„) = ^-« 2 /2r + (,--lW2,.+ l/4 /i(x 2« ); 



fi(z) 



satisfying the equations 
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(4.6) 



/(«) " [«"!]' 

These vertex operators satisfy Eqs. (3.3) and (3.5) and we want to identify them with the operators 
(fri 1 l ' l \u). But they contain a new free parameter uq, whereas the final formulas for correlation func- 
tions like (3.8) and (3.9) must be uo independent. So let us look at the situation more carefully. Any 
u -dependent operator O^' ,l \uo) : Hm — > wt' with Hm = 0„ e 2Z+m+j ^m,n, which is defined in the SOS 
model, is related with the respective uo-independent operator '*) : Hi — > Hy , which is defined in the 
eight-vertex model, by the intertwining property 



where 



O«>*{u )U£\u ) = ug\u )O«' i >, 



Here u is the spectral parameter at the vertical arrows in Fig. 6c, sequences £i,£2,-- - and n, m, ri2, . . . 
stabilize to the respective ground state sequences, and each intertwining vector is divided by its value on 
the ground state configuration, which is denoted by hats. So the traces in both cases must coincide if the 
Mo-dependent operators are properly defined. 
The 'tail operator' A(u) = n n , A(u)£ is 



A ( ) n ™1™2 



= L 



V 










n> 2 








u 


L 






u 


n 


m 






m 


n 2 





with 



n 4 n 3 
m n 2 
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713 ' 



(4.7) 



(4.8) 
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In fact, the r. h. s. of (4.7) is a finite product because all sequences n\, ri2, ■ ■ ■ , rik, ■ ■ ■ stabilize to (— ) k+l by 
definition. Note also that we only need the elements with n' — n e 2Z [or ri\ — e 2Z in Eq. (4.8)]. This is 
because only even number of vertex operators enter into the correlation functions (see Fig. 6c) . So we shall 
imply it from here on. 

The functions L can be easily calculated explicitly: 



V 


n'± 1 


u 




-n')][\(n + n>)] 


n 


nil 




MM 


~n' 


n'T 1 


u 




+ n')][§(n-n')] 


n 


nil 




MM 



By definition 



and hence 



L 



n 


n" 








u 


n 


n> 





A« = 1. 



The counterpart of the operator i s 



for n' - n <E 2Z. 



p«K)= const A( Uo );'[4 



(4.9) 

(4.10) 
(4.11) 

(4.12) 



To prove the defining equations (3.10) we need only prove the properties (3.4-3.5) of the vertex operators. 
The commutation relation (3.3) and the normalization condition (3.5) are proven directly. 236 The most subtle 
point is the proof of Eq. (3.4) for (u, u ) and p^(u ). To do it we need the commutation relation 

between $ e (u; u ) and A(u )- As it shown in Fig. 8, we have the following relation 



V 


s' 








Uq — u 


s 


n 



<&«,A(u )s. 



(4.13) 



Multiplying it by f(u — u — 2)t' £ (u — u )^, taking sum in s, and using Eqs. (2.10), (2.11) we have 

^"^(u; u )A(u ) = A(«o)/(« - uo - 2) & e (u - u - 2)^'$«'. 

, M 



We see that commutation with A shifts u by —2 in each function except &(u)™ and adds the factor [n']/[n]. 
Now pulling <&(«)" through [n] x 4Hn we obtain (3.4). 

It is necessary also to ensure that the spectrum of p^(uo) from (4.12) coincides with the spectrum of 
pW from (3.1). To do it let us calculate the quantity 



4H„ 



-_2Z+m+i 



')= X! MXm,n(a; 4 )- 

ne2Z+m+i 



(4.14) 



As well as we know, M. Jimbo was the first to make this observation. 
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In the Appendix B it is shown that 

x£V) = Kx W (* 4 ) (4-15) 



with 



We expect that traces of operators (4.4) and (4.12) are proportional to [to]' so that this factor cancels in 
the probabilities. This expectation is realized at least in the case of two vertex operators in the trace (see 
Appendix D). In other words, we expect that 



A(u ) Pm = const [m]'x iH ^ 
with ifW(uo) being an operator with the same spectrum as the eight-vertex corner Hamiltonian ffW and 



5. Free Field Representation for the Eight- Vertex Model 

5.1. Free Field Representation for the SOS model 16 

The most efficient way to calculate correlation functions in the SOS model is based on the free field repre- 
sentation. Consider a Heisenberg algebra of operators a n with nonzero integer n and a pair of 'zero mode' 
operators V and Q with the commutation relations 

[P,Q] = -i, [^M = k [k] ^-^ ] ' 6^, with M* = f-^- (5.1) 

The normal ordering operation :. . .: places V to the right of Q and with positive k to the right of a_fc. 
Now introduce the field 

* {z) = -\f^r {Q [V logz) g a i z ~ k - (5 - 2) 

This field enters into the exponential operators 

V(u) = z (r ^ 1)/4r :e iv(z) :, V(u) = z^-V/r^-Mx-^)-^^).^ (5.3) 
and Lukyanov's screening operator 



X(u) = r? _1 e / — V(v 



dv - [v - u + ± - y/2r{r- T)V] 



C u 17r [v - u - \] 

-l = -hl ^ {X 2 ;X^U (x 6 -x\ X ^)oo(x 2r+2 ;x\x 2r )c 
T\ 1[1J X r . 2r~2. rr 2r\ („A.„A ~.2r\ ( , T 2r+4.™4 rr 2r\ 

The contour C„ goes from litoti + j. 

Now the vertex operators are represented as follows: 

:m—n 

- -[^(«). 

:m—n 

- (-) m - n+1 T ^^(«)A( U ). 

[nj 

The corner Hamiltonian is given by 

52 



(5.4) 



(5.5) 



ff = - + gpu^W fl - fcafe - (5 ' 6) 
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The operators (5.5) and (5.6) act in the direct sum of the Fock spaces T m ,n generated by the operators 
with k > from the highest weight vectors \P m ,n) such that 



^fc|-fm,n) (fc > 0), 7~*\Pm,n) ^m,n|-fm,n); Pm,n 



2(r-l) 



r- 1 

2r 



(5.7) 



It is not clear if the spaces (J) n 7i mi „ and (J) n ^-"m.n coincide as representations of the vertex operator algebra. 
At least the numbers of states at each level (according to the H grading) of both spaces Tt m ,n and T m ,n 
coincide. We only assume that traces over both spaces coincide. 

Later we shall need the following fact on the SOS model. The weights (2.6) are invariant under the 
substitution n k — > — rife. So the model is invariant with respect to the change 



n k — » -rift 
m — > —m. 



(5.8) 



The free field representation (5.5) is not invariant under this substitution. So there is another free field 
representation 

•m— n 

= V (u)X(u), 

M (5-9) 
-V(u). 

n . Of course, this representation gives the 



Ti—l / \m— n+1 



^wr 1 = (-) 

In this representation the space Ti m ,n is identified with T- 
same local heights probabilities as the first representation. 



5.2. Free Field Representation for the Operator A(it) 

Now let us turn to the eight-vertex model. The representation of p'M and c^ 1- "-* (u; u ) completely reduces 
to the free field representation of the SOS model by use of Eqs. (4.3) and (4.4). The only unknown object 
is A(ito). To find its representation let us consider the limit uq — > u in the commutation relation (4.13). In 



this limit L 



V 


s' 








u 


s 


n 





L 



-» oo and 

n' n' + l 
nil n 



n 1 n'-l 
n±l n 



u — ► 0. 



We have 



Substitution of (5.5) into (5.10) gives 



V(«)A(« 



,n'-l 



(-) r 



-$( M )™: +1 a( U )^'+ i . 

n' - 1] . 



(5.10) 



\V{u)X(u)A(u) 



n' + l 



[n 1 + 1] 

For n' < n and odd n' — n we have a solution with the initial condition (4.11): 



A« 



n-2fc 



( _ )(m _ n+1)fc [n 2fe] x fc (u) for fc > Q ^ 



(5.11) 



In the Appendix C we prove that this solution satisfies the general commutation relation (4.13). So we may 
consider Eq. (5.11) as the free field representation of A. Now we have to construct representatives for A(u)™ 
with n! > n. However, we know no such representatives in the representation (5.5). To solve this problem 
let us use the second free field representation of the SOS model (5.9), described at the end of Sec. 3. In this 
representation there is a natural representative 



A«+ 2fc = (-)^-" +1 ) fc [n + 2fc] A fc ( M ), for fc>0, 



(5.12) 
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acting from T- m ,-n to ,F_ TOj _ n _2fc. So wc must use the representation (5.5), (5.11) for the contributions 
with n' < n, and the representation (5.9), (5.12) for those with n' > n. 



5.3. Staggered Spontaneous Polarization 

Let us check our construction in the calculation of the staggered polarization in the antiferroelectric regime of 
the eight- vertex model. The staggered polarization is the average 'spin' (e)^ in the state with the condition 
at the infinity i. According to (3.8) it is given by 



(e) W - J2 = Trp (!) M E e Tr (*-e _<) (« " 1. «o) («> u °) P (l) M) • 



Substituting (4.4) we obtain 

;«_(-> 



y— p-(x 2 ; x 4 )^ ^ ^ MX 

' J n 2 £2Z+m+i ni=n 2 ±i 

»0=»1± 1 

Tr^ m ,„ 2 (*(«- 1)-$( U )-A( W0 )^^) 





n 






It — Uq 


m 


"-2 



where 



(-) 



n 3 — m 



A' 



714 


"3 




ni 


«2 





/(«)/(«- i)5>M« -!»(«)£■ 



Applying the trick with the sign change of m and ns we obtain 



( e ) < ' ' — r 1/ ( x ! x )°o E] I n ] ("^n (Tm,n ^-m,n) + + (^m,n ^-m,n)) 
nG2Z+m+i 

n — 1 n — 2 
n — 1 n 



/v.; = a 


n - 
n - 


1 n 
1 n 




- Uq 


m,n 


= Tr^ 




(*(« 


m,n 


= Tr^ 




(*(« 



A-++ = A 



U — Uq 



\7l-l , 



\n-2 4ff„ 



(5.13) 



(5.14) 
(5.15) 

(5.16) 



In Appendix D we perform the direct calculation using the free field representation. The answer is m and 
u independent and is given by the well-known Baxter-Kelland formula 5 



( £ )W =-(-!)« 



f i2L 

edv 4 \ r ' er 

i ^ r ' er J 



{-x 2 - 1 x 2 ) 2 0O {x 2r ;x 2r ) 2 O0 ' 



(5.17) 



This gives an additional 'finite temperature' check of the correspondence of ground states of the eight-vertex 
and SOS models. 



6. Discussion 

We proposed a free field construction for correlation functions of the eight-vertex model. The main scheme 
is the following. We have the bosonic representation for the SOS model, and we know the vertex-face 
correspondence. Using it we can express objects of the eight-vertex models in terms of the objects in the 
SOS model. The only nontrivial point is the representation of the operator A(u) — the unremovable 'tail' 
of the intertwining vectors. The matrix elements A(w)™ _2fe with k > are represented up to a constant 
factor as powers of Lukyanov's screening operator X k (u), and A(u)™ +2k remain undefined. Nevertheless, it 
is possible to calculate necessary correlation functions using the fact that the SOS model admits of two free 
field representations. 
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Our construction is checked in several ways: 

• it supplies solutions to the defining relations for correlation functions; 

• it gives the correct spectrum of the corner Hamiltonian; 

• it gives the correct average staggered polarization. 

Nevertheless there still remain some open problems. The most fundamental problem is why we may 
substitute traces over the configuration space "H m ,n by traces over the Fock T m ^ n l In the case of the six-vertex 
model the identification of the configuration and Fock spaces is based on the quantum affinc symmetry and 
low temperature expansions. 9 ' 25 In the case of SOS model there is no such symmetry algebra, which could be 
defined directly on the lattice. We believe that the solution of this problem could also clarify inconsistency 
of the involutions in the algebra of vertex operators and in the Heisenberg algebra. This inconsistency forces 
us to use two free field representations for vertex operators to find correlation functions. 

An important problem is to prove m and Uq independence of all correlation functions in the framework 
of the free field representation. This would be a strong argument for the consistency of our approach. 

Another group of questions is related with the problem of calculation of form factors. For the six-vertex, 
SOS, and RSOS models we know that there exist type II vertex operators, that describe asymptotic states 
in the models, and allow one to obtain general form factors of local operators. The existence of such vertex 
operators is substantiated by the construction of the elliptic algebra. 26 To find the representation for the 
form factors it seems to be necessary to find a free field representation of generic matrix elements A(u )™ " . 
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Appendix A 

Let us consider a trace like (5.13). The formal series that define this trace is convergent in the region 

|:r~ 4 zo| > \x~ 2 z\ > \z\ > \z \ or Re u — 2 < Rc u — 1 < Re u < Re uq- 

So we must consider the trace for arbitrary values of u and uo as an analytic continuation from this region. 
Besides we know that in the physical region (the regime 77/) 

< u < 1. 

As a consequence of these conditions we have 

-1 < Re(?i - u ) < 0, -2 < Re(-Mo) < 0. 

On the other hand the intertwining vectors are contained in three objects: operator <& £ (u, no) contains 
t £ {u — u ); operator $ e (u — l,tio) contains t*(u — Uo); operator p(uo) contains t e {— u ) and t*(— u ). So 
all the intertwining vectors t e (v)™ , i*(i>)" entering into the correlation function representation are of the 
argument in the region 

-2 < Rev < 0. 

Appendix B 

Here we prove the identity 

£ Mx m ,n(* 4 HM'x w 0A (B.i) 

ne2Z+m+i 
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We have 



ne2Z+m+i 

J2.Pt „ 



E 



nE2Z+m+i 



(x 4 ;x 4 ) 



E 



„4/iL 



ri£2Z+m+i 



(x 4 ;x 4 ) 



-e x2r (x 2n ) 



with 



P 2 1 /n 2 

7 / m.n . -L ( it 

h m,n = — + -A--n 



n 2 <2m+l)n 1 r 9 



We see that the fractional coefficient at n 2 in 2P^ nn is compensated in &h' m n by the fractional coefficient 
that takes its origin in [n]. 

Now let us apply the identity 



- £(-i)V (fc+1 

feez 



We obtain 



rW 



(*) 



{x 4 ;x 4 ) 



kez 

a ,n 2 -(2m+2fe+l) 

n£2Z+m+i 

(fc+m)(fc+m+l) 



r , - ... ^ X 41 21 = x -(k+m)(k+m+l) y^a; 4r 



^2i 



Jez+i(fe+i) 



/ez 



_ ~-(fc+ro)(fc+TO+l) ft j_ J\ _ -(fc+m)(fc+m+l) (g j £ )°° 

Here we used the substitution n — > 2Z + fc + m and then Z — > Z + i for odd fc. The equality in the last line is 
easily obtained by use of the evident identities 



So we have 



(2;p)oo = (z;p )oo(zp;p )«,, (z ;p )«, = (z;p)oo(-z;p) 



(i) _ 1 V^/_ n fe„(r-l)fe(fc+l)-2mfc+^~m _ H' 

(x 2 ;x 4 )-- J ^ 2 - 4 



feez 



(x 2 ;x 4 ) c 



q. c. d. 

Note that this proof, as well as the consideration in Appendix D, is valid for to, n e Z + 5 with any (5. 



Appendix C 

In this Appendix we prove that the operators (5.11) satisfy the commutation relation (4.13). Namely 
consider the quantity 



d( !1 , M „)=A«r« + -i 

- L 



n+l-2fc n+2-2fc 
n + 1 n 



n + l-2k n-2k 
n+1 n 



uq — u 



it — u 

jl+l-2/c a /„, \n- 2fe 



n+l-2fc A / , N,n+2-2fc 



(C.l) 



1G 



We have to prove that d(u, uq) = from the free field representation inputs. Substituting Eq. (5.5) we obtain 

d(u, «o) ~ [n + 1 - 2k]X k {u )V(u) + (-l) fc [U ° - - r + n + ) - m V{u)X{u)X k -\u ) 

[u - u\ 

fe [n -^ + fc][n+l-fc] fc 
[« - uj 

Here the sign ~ only means that two functions differ by a factor which is allowed to be any nonzero function 
of all variables. Substituting Eq. (5.4) and taking into account that 



[ui - u 2 + |] 

[«i -u 2 -h] 



V( Ul )V(u 2 ) = 1 ^V(u 2 )V( Ul ) (C.2) 

we obtain 



d(u,u )~ J k V(u)V(v k ) . . . V(vi)f(u, u ;vi,..., v k ) 



d k v 
(2m) k 
with 

t, \ a fr , 1 on A K ~ u + 3 - n - 1 + 2 (i - 1)] K - w + \\ 
f(u, uo; «i, . . . , ffc) = Sym<^ [n + 1 - 2 A: I I n - ft 

[u - u + n + 1 - k][k] [vk - u - n + 2k - §] tj- 1 [w, - u + \ - n + 2(i - 1)] 
[tio -u] [v k -u - §] [„. - Uo - I] 

[up - it + fc] [n + 1 - k] yr [vj - u + \ - n + 2(i - 1)] | 



[tlO - tl] [Vj - «0 - I 

Here the operation Sym is defined as 

SymF(ui, ... ,v fe ) = ^ F (v a{1) ,... ,v a{k) ) J| /i(tv w - t; CT(j) ) 

with /i(u) being defined as 

V( Vl )V{v 2 ) = h( Vl - v 2 )V{v 2 )V{v 1 ), h(v) = ^— ^ 

[v + 1] 

This 'symmetrization' of the function / corresponds to symmetrization of the integrand. 
So we have to prove that f(u, u a ; v\, . . . , v k ) = 0. 
To calculate / let us use the identity 27 

i— 1 L J i<j 1 J 1 i—1 

with [fc]! = Yli=i W- We apply this function directly to the first and last terms, but in the second term we 
only apply it to fc — 1 variables, and symmetrize in the remaining variable by hand. We have 

f(u,u ;vi, . . .,v N ) ~ g(u) 

with g(u) being a doubly periodic function of u (and all other variables): 

, . _ [n + 1 - 2fc] [u - it] yr[v i -u -n + k-\}[v i -u+ \] 
3 U [n-k + 1] [uo - u + k]*^[vi - uo - n + k - \][vi - u - \] 

[uq - u + n+ 1 - fc] [1] y-r [vi - u - n + k - §] 



x 



[tio - u + fc] [n - k + 1] [vi - uo - n + k - |] 

y, / [vj - u - n + 2fc - |][^ - up - |] „ [vj - Vj + 1 ] 
£t I ^ - « - 5]^ - uo - n + fc - |] 1A K - «j] 
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It is easy to check that 

g(uo + n + 1 — k) = 0. 

So it is enough to prove that the function g(uo) is constant. It has poles at the points (up to periods) 

u = uq + k, 
u = vi- 1/2. 

It is straightforward to check that the residue at u = vi — | vanishes. For the residue at u = u + k we have 
Res g H^ g0 M = l-^n 5^' M °'fi| vf r hz^^ T\h^i±}\ 

u=u 0+ k yy [k] 11 [ Vi - U - k + i] ^ ^ [ Vj - U - k - i] 11 [«i - 

It is easy to check that go(vi — k — |) = 0. For the residues we have 



, 50 WO (Wj(l) = 1- 77 TT V M -T-Tj j u TT 



To proceed by induction introduce the functions 

We suppose that . . . , is an arbitrary permutation of the numbers 1, . . . , A; (i. e. I € Sk)- We have 

9s{vi(s+i)) = 0, Res o a (wi( s )) ~ 5 s+ i(wi( 5+ i)), 9k{vi{k)) = 0. 

" 1(s) =i> i(s+1) +i 

Hence it follows that all g s and / vanish. This completes the proof in the case s = n + 1. The proof in the 
case s = n — 1 is similar though more cumbersome. 



Appendix D 

Here we show that the free field representation for the eight-vertex model results in the Baxter-Kelland 
formula for the staggered polarization. 

As the initial point of the calculation we take Eq. (5.16). For the coefficients K we have 



K 



K 



V 


n'±l 




h 4 (u i 




n 


nil 




hi(u) 


V 


n'T 1 


u 


/l 4 (w i 


^f)h 4 (J^f) 


n 


nil 




hi(u) 



for n' — n e 2Z, with 
Hence 



7T l 

e 4&r 6 
er 



W 17T 

r ' er 



hi(u) = [u], h 4 (u) = ie I7r " [u - i7r/2e]. 

+ _ _ /t 4 (n - l)fe 4 (^) 
++ _ h 4 (n-l + u)h 4 (0) 



(D.l) 



(D.2) 
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The traces T+ and T++ can be expressed as follows 

T " l - n = ~[n-l][n] Tm < ni * U,U) ' T ^ n = [n-l][n] Tm - n ^ Uo) 

with 

T m , n (u,u ) = Trr m jV(u-l)V(u)X(u )x AH ) 

= e [ ^-H m n (u,u ;v)G(v - u). (D.3) 

J 171" 

Here G(v) is the contribution into the trace in the second line from the oscillators a n : 



= i 7 ^Le«^-«^(--i) a ^4| I C = ^, (D.4a) 



G( V ) 



r y - 1 AA=ix- I ^(x 2 ; a ; 2 ) 2 0O (x 2r ;x 2 '') 0O , ^(u) = ^ e^u; iir/e). (D.4b) 

The factor A/" originates from the traces of individual exponential operators, whereas the rest originates from 
traces of pair products in the generalized Wick theorem. 

The factor H min (u,Uo;v) denotes the contribution of the 'zero mode': 



(D.5) 



flm,n(u,uo;«) = ^P ro ,„ * e-'V— y ( x - 2 z)-V — e -'V— y z -V — 

[«-«o-l/2] 

_ ^_^„-m e -2£^i(t)-«+i)-2£(ti-M+i)m-£ 7 ^ T m 2 -e^in 2 +2emn+2e(t;-ti+i)r 

/ii (u - u + i - n) 
hi(v - u - \) 

It is convenient to calculate the quantity 

H${u,u ;v)= T^^{K+ + H m . n {u,u ]v) - K+^H m . n {u,u]v)) 

— 26-^—5- (v— u+j) — 2e(v— u+^)m — e—^jm 2 - e — n 2 ~2emn+2e — (u-u+|)n 

^ hi(u - u ) hx{n - 1) hi(v - u - 5 ) - w - 5) 

x {ft-i(f — tt — i) /ii(i> — Mo + \ — n) hi{n — 1 + u — uq) /i 4 (0) 

— — i*o — 5) h\(v — u + i — n) /i 4 (n — 1) /i 4 (w — u )}. (D.6) 

Substituting the identity 

4 (2ar)0 4 (2j/)0i(2z)0i(2i) 

= 6*4(2; + y + z + t)6i{x + y — z — t)0\(x — y — z + t)0\{x — y + z — t) 

+ 6 4 (x + y + z- t)6i{x + y- z + i)6>i(x - y + z + t)6i{-x + y + z + t), (D.7) 

we obtain 

H$ (U, U ; v) = (-)* e -2e^i(i,-u+i)-2£(«-«+i)m-£ 7 ^ T m 2 - e ^in 2 +2em«+2e^i(t,-«+i)n 
n 

/i 4 (v - uq — 5) /i 4 (w — w + 5 - n) 
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(D.8) 



From the summand we can extract an n-dependent factor 

\ r er J \j tt \ n ir / 

with w = v — u + i. Expanding the last theta function in a series, we obtain 

^(i) _ €r c -ers 2 -2ews~4ep 2 +4:emp+4ewp+4esp-^w 2 
ne2Z+m+i n p£ z+i(m+i) s£Z+i 

Now we can split the sum into two terms: the first with even s — m — \ and the second with even s — m+ 5. 
Each term factors after change of summation variables (p, s) to (k, I) so that s = m ± \ + 21 and p = k + I. 
Each sum in k and I gives a theta function. 

H£Hu,Uo;v) = hYe hr-2^ {v -u +i)+ ^ {v -u +i f h 4 (v - no |) 

— uo — |J — u — 5) 

„ (™ - ^) , ^) (if („ + i) , if) 

Here the subscript 3 or 2 depends on evenness of m + i. Making the module conjugation transformation we 
get 



H$(u,u ;v) = — ^=< -< 



c ier-2e^(v-u+i)+e^l (v -u+^ h 4 { V - Up - ±) 

4\/r - 1 - wo - \) h\(v — u — |) 



fir 1 yk \ „ ( w 1 in 
x < f 3 m ; — ; ] 63 4 h -; — 

' ^V 2 ^- 1 4'4e(r-l)y ' \2 4' 4e 

+ 03 (2 ~ m + i ; i^T) ) ^ ( 2 - I ; 47) }• (D ' 9) 

Here the subscript 3 of the second and forth theta functions corresponds to even m + i and the index 4 to 
odd one. In principle, the above resummation procedure is the same as that of the Appendix B. 
Calculating Hm — H^ m it is convenient to use the identity 

6 3 (u + ±; t) - 3 (u - \- t) = -2 9 1 (2u; At). 

Applying it twice, we get 



h^v-uo- \)hi_{v-u- |) 



(D.10) 



The factor [m]' here will cancel the respective quantity in the denominator in Eq. (5.16). The u depending 
theta functions here will cancel theta functions in G(v). 
Gathering all results we obtain 

(e)W=i( _ ir (^!k £ J ^ { H^(u,uo;v)-H^ m (u,uo;v))G(v-u) 
(x 2 -x 2 ) 2 

= -(-l) V 2 2T2 ^ C - 11 ) 



/oo 
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where 

I=e *L 6 <—>~) (D . 12) 



Co 



17T 



\ r ' er y 



To calculate the integral note that 

0i(u+l)=-0i(u), 9 A (u + l) = 9 4 (u), 

and 



Co i?r fli f ^+ r ~ 1 / 2 . i^'l "7 Cr i7T g ( v = l/2 .m 

ii r i er J 1 V r ' er 

But the sum of (D.13) and (D.12) is an integral over a closed path enclosing the pole v — 1/2. Hence 



Using the multiplicative form of the theta functions we easily obtain 
Substituting (D.14) into (D.ll) we obtain the Baxter-Kelland formula (5.17). 



(D.14) 
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